Int. J. Mol. Sci. 2014, 15, 7987-8003; doi:10.3390/ijmsl5057987 



OPEN ACCESS 



International Journal of 

Molecular Sciences 

ISSN 1422-0067 

www.mdpi.com/joumal/ijms 

Article 

Endothelial Nitric Oxide Synthase (eNOS) 4b/a 
Gene Polymorphisms and Coronary Artery Disease: 
Evidence from a Meta-Analysis 

Yujiao Yang ^, Kang Du ^, Zhengxia Liu and Xiang Lu * 

Department of Geriatrics, the Second Affiliated Hospital, Nanjing Medical University, 
Nanjing 210029, Jiangsu, China; E-Mails: yangyujiao626@126.com (Y.Y.); 
dukang@njmu.edu.cn (K.D.); liuzx0503@gmail.com (Z.L.) 

^ These authors contributed equally to this work. 

* Author to whom correspondence should be addressed; E-Mail: luxiang66@njmu.edu.cn; 
Tel.: +86-25-5880-2883; Fax: +86-25-5850-9994. 

Received: 25 March 2014; in revised form: 24 April 2014 / Accepted: 29 April 2014 / 
Published: 7 May 2014 



Abstract: A variety of studies have suggested that the 4b/a polymorphism in the endothelial 
nitric oxide synthase {eNOS) was associated with coronary artery disease (CAD) risk. 
However, the data remain conflicting. The aim of the present meta-analysis was to estimate the 
overall association between risk of CAD and eNOS 4b/a polymorphism. Case-control, cohort 
or cross-sectional studies evaluating the association between eNOS 4b/a polymorphism and 
CAD susceptibility were systematically identified in PubMed up to 31 October 2013. 
Pooled odds ratios (OR) and corresponding 95% confidence intervals (CIs) were calculated 
to assess the association in overall and subgroup analyses. A total of 10,617 cases and 
8302 controls from 37 studies were included in the study. The results of overall analysis 
revealed significant positive associations between CAD risk and eNOS 4b/a polymorphism 
in homozygote comparisons (OR = 1.47, 95% CI = 1.16-1.87), heterozygote comparisons 
(0R= 1.14, 95% CI = 1.02-1.27) and dominant models (OR = 1.18, 95% CI = 1.06-1.33). 
In subgroup analyses, similar associations were identified in African individuals, 
as determined using population-based source subgroups and noted in small-and-moderate 
sample size subgroups (case sample size or control sample size <500). The current 
meta-analysis revealed that eNOS 4b/a polymorphisms could be a risk factor for developing 
CAD, particularly in African populations and population-based subgroups. 
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1. Introduction 

Coronary artery disease (CAD) also known as coronary heart disease (CHD), is the leading cause of 
death and disability worldwide [1,2]. CAD is associated with genetic and environmental factors, as well 
as their interactions [3]. Confirmed risk factors for CAD include hypercholesterolemia, hypertension, 
smoking, and diabetes [4]. However, in addition to modifiable risk factors, genetic factors can also 
predispose individuals to CAD. It has been estimated that genetic risk factors explain approximately 
20%-60% of CAD cases [5]. 

Nitric oxide (NO) is an important atheroprotective mediator that helps mediate endothelium-dependent 
vasodilatation. Abnormalities in NO generation could play an important role in the pathophysiology of 
CAD [6]. NO is synthesized from L-arginine by the action of nitric oxide synthase (NOS). There are at 
least three isoenzymes of NOS: inducible NOS, neuronal NOS, and endothelial NOS (eNOS) [7]. 
The eNOS gene is located on chromosome 7q35-q36 and is sized 21 kilobases (kb) [8]. A number of 
studies have suggested that polymorphisms in the eNOS gene affected NO availability and were associated 
with CAD morbidity. One of the most clinically relevant eNOS variants forms 27-basepair (bp) tandem 
repeats in intron 4 (4b/a). The eNOS 4b allele contains five repeats of the 27 bp sequence, whereas 
eNOS 4a, which is the rare allele, carries a deletion of one of the first three of these repeats [9]. 

The association between eNOS 4b/a gene polymorphisms and CAD remains controversial, and the 
results of two previous meta-analyses were inconsistent [10,11]. In 2004, Casas et al. performed a 
meta-analysis of 16 studies and found an increased risk of ischemic heart disease (IHD) in individuals 
carrying the homozygous a allele compared with the b allele [10]. However, a subsequent meta-analysis 
in 2010 reported conflicting data, where eNOS 4b/a did not significantly increase the risk of developing 
CAD [11]. A large number of new case-control studies assessing the association between eNOS 4b/ a 
polymorphisms and the CAD risk have been published since 2010. Therefore, we performed a 
meta-analysis using a more complete database to clarify the association between eNOS 4b/a 
polymorphisms and CAD risk. 

2. Results and Discussion 

2.1. Characteristics of Included Studies 

A total of 373 papers were initially identified during the literature search and from other sources; 
37 of these met the inclusion criteria for this analysis [9,12-47]. A flow chart describing the literature 
search and study identification is shown in Figure 1. The meta-analysis of the 4b/a polymorphism 
included 10,617 CAD cases and 8302 controls. The study characteristics are summarized in Table 1. The 
number of subjects in the studies varied considerably (range from 40 to 1265 in the case group and 
from 34 to 620 in the control group). Among the 37 case-control studies, 1 1 studies included Asian 
population [15-25], and 26 assessed non-Asians [9,12-14,26-47]. Controls were mainly matched in 
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terms of gender and age. Fourteen studies were hospital-based [16-18,20,21,25,28,37^0,42,43,45], and 
20 were population-based [9,12-15,19,22,23,26,27,29-34,36,41,46,47]. Thirty studies were performed 
using <500 subjects in the case group [12-27,29,33-36,38^7] and 33 were conducted with <500 
subjects in the control group [9,12-15,17-29,32,34^7]. Furthermore, there were 20 studies 
of CAD [9,18,20,22,23,25,29,30,32,34^0,42^4,46], four of CHD [17,26-28], three of acute 
myocardial infarction [15,19,45], six of myocardial infarction (MI) [12-14,16,31,41], and four of acute 
coronary syndrome [21,24,33,47]. 

Figure 1. Flow diagram of study identification. 
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2.2. Meta-Analysis Results 

The relationship between eNOS 4b/a polymorphisms and CAD risk are summarized in Table 2. 
Overall, significant positive associations between CAD risk and eNOS 4b/a polymorphism were identified 
in the homozygote comparison {aa vs. bb, OR = 1.47, 95% CI = 1.16-1.87, /»heterogeneity = 0.085), the 
heterozygote comparison (ab vs. bb, OR = 1.14, 95% CI = 1.02-1.27, /^heterogeneity < 0.001) and the 
dominant model (ab + aa vs. bb, OR=\.\ 8, 95% CI = 1 .06-1 .33, /^heterogeneity < 0.001). 
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Table 1, Study Characteristics of genotypes in coronary artery disease (CAD) cases and controls in the analysis of endothelial nitric oxide 
synthase (eNOS) 4b/a polymorphism. 

Genotype distribution 



Author 


Year 


Country 


Ethnicity 


aource oi 
controls 


End points 






Control 






Case 














Control 


Case 


aa 


ab 


bb 


aa 


ab 


bb 


Wang et al. 


1996 


Australia 


Caucasian 


PB 


CAD 


153 


549 


1 


50 


102 


18 


122 


409 


Hibi et al. 


1998 


Japan 


Asian 


PB 


AMI 


357 


226 


5 


68 


284 


4 


48 


174 


Ichihara et al. 


1998 


Japan 


Asian 


HB 


Ml 


550 


455 


7 


97 


446 


9 


107 


339 


Odawara et al. 


1998 


Japan 


Asian 


HB 


CHD 


122 


42 


0 


19 


103 


0 


11 


31 


Hooper et al. 


1999 


American 


African 


PB 


MI 


185 


110 


15 


68 


102 


14 


52 


44 


Nakagami et al. 


1999 


Japan 


Asian 


HB 


CAD 


34 


40 


1 


9 


24 


1 


18 


21 


Park et al. 


2000 


Korea 


Asian 


PB 


AMI 


206 


121 


1 


49 


156 


5 


23 


93 


Fowkes et al. 


2000 


UK 


Caucasian 


PB 


CHD 


300 


137 


5 


64 


231 


4 


37 


96 


Pulkkinen et al. 


2000 


Finland 


Caucasian 


PB 


CHD 


110 


308 


5 


26 


79 


9 


96 


203 


Sigusch et al. 


2000 


Germany 


Caucasian 


HB 


CHD 


413 


625 


11 


115 


287 


12 


159 


454 


Granath et al. 


2001 


Australia 


Caucasian 


PB 


CAD 


620 


567 


14 


158 


448 


11 


138 


418 


Alvarez et al. 


2001 


Spain 


Caucasian 


PB 


CAD 


300 


170 


6 


72 


222 


2 


37 


131 


Hwang et al. 


zuuz 


Taiwan 


Asian 


rlrJ 




70 


149 


A 
U 


1 A 


JO 


1 

i 


jL 


1 1 


Gardemann et al. 


imi 


Germany 


Caucasian 


PB 


MI 


528 


1265 


13 


144 


371 


25 


306 


934 


Cine et al. 


2002 


Turkey 


other 


PB 


MI 


306 


207 


2 


55 


249 


9 


55 


143 


Park et al. 


2004 


Korea 


Asian 


HB 


ACS 


142 


164 


7 


30 


105 


4 


21 


139 


Agema et al. 


2004 


Netherlands 


Caucasian 


PB 


CAD 


466 


752 


8 


11 


381 


12 


195 


545 


Letonja et al. 


2004 


Slovenia 


Caucasian 


PB 


CAD 


109 


151 


6 


30 


73 


6 


41 


104 


Fatini et al. 


2004 


Italy 


Caucasian 


PB 


ACS 


537 


477 


14 


138 


385 


24 


138 


315 


Milutinovic et al 


2005 


Slovenia 


Caucasian 


NM 


CAD 


188 


215 


10 


58 


120 


11 


60 


144 


Rao et al. 


2005 


America 


Caucasia + 
African 


HB 


CAD 


50 


144 


1 


17 


32 


6 


40 


98 


Matyar et al. 


2005 


Turkey 


other 


HB 


CAD 


133 


133 


1 


35 


97 


5 


47 


81 


Salimi et al. 


2006 


Iran 


other 


NM 


CAD 


158 


141 


1 


29 


128 


3 


41 


97 
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Table 1. Cont. 











Source of 
controls 




Sample size 




Genotype distribution 




Author 


Year 


Country 


Ethnicity 


End points 




Control 






Case 














Control 


Case 


aa 


ab 


bb 


aa 


ab 


bb 


Agirbasli et al. 


2006 


Turkey 


other 


HB 


CAD 


100 


100 


0 


21 


79 


4 


21 


75 


Kim et al. 


2007 


Korea 


Asian 


PB 


CAD 


211 


147 


0 


40 


171 


2 


35 


110 


Kios al. 


2007 


Brazil 


Caucasia + 
African 


T T r) 

HB 


CAD 


ZD / 




1 o 

12 






1 8 
1 o 


1 01 

lUl 


204 


Sampaio et al. 


2007 


Brazil 


other 


HB 


AMI 


103 


115 


4 


32 


67 


4 


29 


82 


Vasilakou et al. 


2008 


Greece 


Caucasian 


PB 


CAD 


161 


209 


5 


39 


117 


9 


60 


140 


Meluzin et al. 


2009 


Czech rep 


Caucasian 


HB 


CAD 


89 


321 


0 


26 


63 


10 


103 


208 


Kind et al. 


2009 


Czech rep 


Caucasian 


HB 


CAD 


222 


939 


9 


63 


150 


35 


272 


632 


Bae et al. 


2010 


Korea 


Asian 


PB 


CAD 


196 


192 


0 


35 


161 


2 


42 


148 


Gumrajan et al. 


2010 


India 


Asian 


NM 


ACS 


100 


106 


4 


19 


77 


12 


36 


58 


Jaramillo et al 


2010 


Chile 


other 


PB 


CAD 


112 


112 


1 


16 


95 


3 


12 


97 


Jemaaa et al. 


2012 


Tunisia 


African 


PB 


MI 


250 


310 


4 


61 


185 


15 


105 


190 


Da Costa Escobar Piccoli et al. 


2012 


Brazil 


other 


PB 


ACS 


108 


132 


6 


34 


68 


14 


31 


87 


Kallel et al. 


2013 


Tunisia 


African 


PB 


MI 


225 


303 


3 


58 


164 


15 


101 


187 


Name et al. 


2013 


India 


Asian 


HB 


CAD 


121 


160 


2 


42 


77 


5 


50 


105 



Abbreviations: eNOS, endothelial nitric oxide synthase; HB, hospital-based; PB, population-based; NM, not mentioned; 
CAD, coronary artery disease; CHD, coronary heart disease; MI, myocardial infarction; AMI, acute myocardial infarction; and 



OR, odds ratios; CI, confidence intervals; 
ACS, acute coronary syndrome. 



Table 2. Pooled odds ratios (OR) and 95% confidence intervals (CI) of the association between eNOS 4b/a polymorphism and CAD. 



Total and 
subgroups 



Studies 



Homozygote comparison 
{aa vs. bb) 



Heterozygote comparison 
{ab vs. bb) 



Dominant model 
{aa + ab vs. bb) 



Recessive model 
{ab + bb vs. aa) 



OR 

(95% CI) 



I' (%) 



OR 

(95% CI) 



I' (%) 



OR 

(95% CI) 



f (%) 



OR 

(95% CI) 



I' (%) 



Total 
MI 



37 
6 



1.47 
(1.16-1.87) 

2.30 
(1.17^.51) 



0.085 
0.020 



25.5 
62.6 



1.14 
(1.02-1.27) 

1.42 
(1.07-1.89) 



<0.001 
0.002 



56.7 
74.2 



1.18 
(1.06-1.33) 

1.50 
(1.09-2.07) 



<0.001 
<0.001 



61.1 
80.7 



0.72 
(0.58-0.89) 

0.56 
(0.39-0.82) 



0.213 15.4 
0.054 54.1 
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Table 2, Cont. 



Total and 
subgroups 



Studies 



Homozygote comparison 
(gg vs. bb) 



Heterozygote comparison 
(ab vs. bb) 



Dominant model 
(gg + ab vs. bb) 



Recessive model 
(ab + bb vs. aa) 



OR 

(95% CI) 



I' (%) 



OR 

(95% CI) 



I' (%) 



OR 

(95% CI) 



OR 

"^^"-^"> (95% CI) ^ "^^"-^°> 



Ethnicity 



Caucasian 
Asian 
African 
Other 



16 
11 
5 
7 



1.03 
(0.80-1.31) 

1.82 
(1.03-3.20) 

2.93 
(1.71-5.03) 

2.65 
(1.37-5.13) 



0.621 
0.316 
0.821 
0.369 



0.0 
13.8 
0.0 
7.8 



1.04 
(0.90-1.20) 

1.21 
(0.96-1.54) 

1.46 
(1.16-1.83) 

1.17 
(0.84-1.62) 



0.006 
0.027 
0.353 
0.033 



53.3 
50.5 
9.4 
56.4 



1.04 
(0.91-1.19) 

1.26 
(0.99-1.61) 

1.58 
(1.27-1.96) 

1.29 
(0.94-1.77) 



0.007 
0.011 
0.370 
0.033 



53.0 
56.3 
6.5 
56.3 



0.98 
(0.77-1.25) 

0.60 
(0.36-0.99) 

0.41 
(0.24-0.70) 

0.39 
(0.21-0.71) 



0.679 



0.434 



0.684 



0.511 



0.0 
0.4 
0.0 
0.0 



Source of control 



HB 
PB 

NM 



14 
20 
3 



1.11 

(0.78-1.56) 

1.64 
(1.18-2.27) 

2.05 
(0.66-6.32) 



0.509 
0.070 
0.111 



0.0 
33.9 
54.4 



1.04 
(0.88-1.23) 

1.16 
(1.00-1.35) 

1.54 
(0.80-2.98) 



0.072 



<0.001 



0.011 



38.2 
61.8 
77.6 



1.07 

(0.89-1.27) 

1.21 
(1.05-1.40) 

1.63 
(0.80-3.30) 



0.027 



<0.001 



0.003 



46.8 
63.4 
82.6 



0.90 

(0.64-1.26) 

0.65 
(0.48-0.89) 

0.58 
(0.24-1.40) 



0.613 0.0 
0.113 28.7 
0.232 31.6 



Case size 



>500 

<500 



6 

31 



0.88 
(0.62-1.24) 

1.75 
(1.35-2.28) 



0.666 

0.262 



0.0 

13.2 



0.96 
(0.74-1.25) 

1.20 
(1.07-1.35) 



<0.001 
0.011 



78.7 
40.7 



0.97 
(0.76-1.22) 

1.26 
(1.11-1.42) 



0.001 

0.002 



75.9 
48.0 



1.14 
(0.80-1.60) 

0.61 
(0.48-0.77) 



0.642 0.0 

0.453 0.8 



Control size 



>500 
<500 



' (0.72-2.05) (0.85-7.36) 



33 



1.56 
(1.19-2.05) 



0.113 23.9 



1.15 
(1.01-1.31) 



<0.001 55.6 



1.09 
(0.84-1.42) 

1.20 
(1.06-1.36) 



0.007 
<0.001 



75.1 
59.0 



0.84 
(0.52-1.34) 

0.68 
(0.53-0.88) 



0.214 33.1 
0.233 14.8 



Abbreviations: eNOS, endothelial nitric oxide synthase; OR, odds ratios; CI, confidence intervals; Ml, myocardial infarction; HB, hospital-based; PB, population-based; 
and NM, not mentioned. 
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In a stratified analysis by the source of the control population, significantly increased CAD risk 
was detected in the population-based subgroup of all the three genetic models (OR = 1.64, 
95% CI = 1.18-2.27,^3 = 0.070 for the homozygote comparison; OR = 1.16, 95% CI = 1.00-1.35, 
p < 0.001 for the heterozygote comparison; and OR = 1.21, 95% CI = 1.05-1.40, p < 0.001 for the 
dominant model). There was no significantly increased CAD risk in the hospital-based subgroup 
(OR = 1.11, 95% CI =0.78-1.56, p = 0.509 for the homozygote comparison; OR = 1.04, 
95% CI = 0.88-1.23,;? = 0.072 for the heterozygote comparison; and OR = 1.07, 95% CI = 0.89-1.27, 
p = 0.027 in the dominant model) (Figure 2). 

Figure 2. Forest plot of CAD risk associated with the eNOS 4b/a polymorphism by source 
of control (homozygote comparison). 



Study 
ID ' 



OR (95% CI) 



population-based 
Wangetal. (1996) 
Hibi et al. (1998J 
Hooper et al (1999) 
Park et a! (2000) 
Fowkes et al (2000j 
Pulkkinen et at [20()0) 
Granath et a! (2001) 

Alvarez et al (200 lj — 
Gardemann et al (2002) 
Cine eta! (2002) 
Agema et al (2004) 
Lipoma et al 12004) 
Fatini et al. (2004) 
Kim et al. (2007) 
Vasilakou et al. (2008) 
Bae et al. (2010y 
Jaramillo et al 12010) 
Jemaaa et al (2012) 
da Costa Escobar F^iccoli et al. (2012) 
Kallel et al. (2013) 

Subtotal (l-squared = 33.9%, p = 0.070) 

hospital-based 
Ichihara et al. (1998) 

Nakaqami et al. (.1999) 

Sigusch et a! (2000) 

Hwang. et al. (2002). 

Park KW et al (2004) — 
Rao et al. (2005) 
Matyar et al. (2(J05) 
Agirbasli et al. ^0(16) 
Rios et al (2007) 
Sampaio et al. (2007) 
Meluzin et al. (2009 
Kind et al (20D9) 
Name et al (2013) 
Odawara et al (1998) 
Subtotal (l-squared = 0.0%, p = 0.509) 

not mentioned 
Milutinovic et al. (2005) 
Salimi et al (2006) 
Gururajan et al {2(J10) 
Subtotal (l-squared = 54.4%. p = 0.111) 

Overall (l-squared = 25 5%. p = 0 085) 
NOTE Weights are from rpndom effects analysis 



7 
1. 
0 
2 
7 
1. 
5 




4.49 
1.31 
2.16 
8.39 
1.92 
0 70 
0 84 
0 56 
0 76 
"84 
.05 
70 
.10 
76 
.50 
44 
2 94 
3.65 
1.82 
4.39 
1.64 



1.69 
1.14 

0 69 

1 45 
043 
1.96 
5.99 
9.48 
1.21 
0 82 
6.40 
0 92 
1.83 



0 59. 
0 35, 
0 96, 
,0.96, 

ro.5i, 

^0 23, 
0 38, 
Oil, 
0.39, 
1,67, 
0.42, 
0.22, 
1.07. 
0.37. 
049, 
0.26, 

0 30, 
1.19, 
0.67, 
1.25, 

1 18. 



(0.62, 
^0 07, 
0 30, 
0 06, 
0 12, 
0.23, 
0,69, 
,0.50, 
ro 57, 
ro.20, 
ro.37, 

^0.43, 
0.35, 



34 02) 
4 93) 
4 86i 
72.90) 
7.32 




114 

28 75)' 
11.21 
5.00) 
15.42) 

2 27) 

Am 

19.42) 
1 58) 
36 28) 

1 5U 
16 89 
52 30 
179.02) 
2.59) 

3 39l 
110.67) 
1.96) 
9.70 



11 (0 78: 1.56) 



0.92 (0 38, 
3.96 r0.41, 
3.98 ri.22, 
2.05 0.66, 



2.23) 
38 65) 
12.98 
6.32) 



1.47 (1.16, 1.87) 



0.00559 



179 



The data were also stratified into Asian and non-Asian subpopulations. There was a statistically 
significant association in the non- Asian subgroup in the homozygote comparison (OR = 1.41, 
95% CI = 1.08-1.84, j3 = 0.080) and the dominant model (OR = 1.16, 95% CI = 1.02-1.32, < 0.001). 
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Non-Asians were then further stratified into Caucasians, Africans, and others. The populations of an 
American study and a Brazilian study were both composed of two ethnicities: Caucasians and 
Africans [39,40]. The different ethnic descents in these studies were considered to be two independent 
studies during sub-sfratified analyses. Statistically significant associations were noted in African 
populations in three genetic models (homozygote comparison OR = 2.93, 95% CI = 1.71-5.03, 
/» = 0.821; heterozygote comparison OR = 1.46, 95% CI= 1.16-1.83,/? = 0.353; and the dominant model 
OR = 1.58, 95% CI = 1.27-1.96,/? = 0.370). Interestingly, significant associations were also found in 
Asian (OR = 1.82, 95% CI = 1.03-3.20, p = 0.316) and "other" (OR = 2.65, 95% CI = 1.37-5.13, 
p = 0.369) populations during the homozygote comparison. In confrast, no statistically significant 
associations were identified in Caucasian populations (OR= 1.03, 95% CI = 0.80-1.31,/? = 0.621 in the 
homozygote comparison; OR = 1.04, 95% CI = 0.90-1.20,/? = 0.006 in the heterozygote comparison; 
OR = 1.04, 95% CI = 0.91-1.19,/? = 0.007 in the dominant model) (Supplementary Figure SI). 

The data were then also stratified according to case sample size into large (>500 total cases) and 
small-and-moderate case sample (<500 total cases) subgroups. Significant associations between 
eNOS 4b/a polymorphisms and CAD risk were observed in the small-and-moderate case sample 
subgroup (OR = 1.75, 95% CI = 1.35-2.28, p = 0.262 in the homozygote comparison; OR = 1.20, 
95% CI = 1.07-1.35,/? = 0.011 in the heterozygote comparison; and OR = 1.26, 95% CI = 1.11-1.42, 
p = 0.002 in the dominant model) but not in the large case sample subgroup. The data were then 
stratified into large control (>500 control individuals) and small-and-moderate control sample 
(<500 confrol individuals) subgroups. Similar to the data from the case sample size analyses, significant 
associations were observed in the small-and-moderate confrol sample subgroup (OR = 1.56, 
95% CI = 1.19-2.05,/? = 0.113 for the homozygote comparison; OR = 1.15, 95% CI = 1.01-1.31, 
p < 0.001 for the heterozygote comparison; and OR = 1.20, 95% CI = 1.06-1.36, p < 0.001 in the 
dominant model), but not the large case sample subgroup. 

Finally, a subgroup analysis for MI was performed. There were significant associations between 
eNOS 4b/a poljmiorphisms and CAD risk with the homozygote comparison (OR = 2.30, 
95% CI = 1.17^.51,/? = 0.020), heterozygote comparison (OR = 1.42, 95% CI = 1.07-1.89,/? = 0.002), 
and dominant model (OR = 1.50, 95% CI = 1.09-2.07,/? < 0.001) in patients with MI. 

2.3. Evaluation of Heterogeneity 

There was minor heterogeneity among studies in the recessive model {ab + bb vs. aa, 
/^heterogeneity = 0.213, / = 15.4%)), while the results of heterogeneity tests were significant in the other 

2 2 

three genetic models {aa vs. 66,/?heterogeneity = 0.085, / = 25.5%; ab vs 6Z',/?heterogeneity < 0.001, / = 56.1%; 
ab + aa vs. 66, /?heterogeneity < 0.001, / = 61.1%). The random-effects models were used in all the four 
genetic models. Meta-regression analyses were then performed to evaluate the extent to which different 
variables explained the heterogeneity. The results revealed that the heterogeneity could be explained by 
case sample size {p = 0.015) and total sample size {p = 0.015) in the homozygote comparison. However, 
population ethnicity, year of publication, source of control population, total sample size, and outcome 
were not statistically correlated with heterogeneity (/? > 0.05). 
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2.4. Sensitivity Analysis and Publication Bias 

Sensitivity analysis was performed by omitting one study at a time. No significant differences in the 
data were observed, indicating that the results were statistically rehable (Supplementary Figure S2). 
The publication bias of the studies was assessed by performing the Begg's funnel plot and Egger's test. 
There was no obvious publication bias in the heterozygote comparison (ab vs. bb, Begg's test p = 0.754, 
Egger's test p = 0.473), and the dominant model {ab + aa vs. bb, Begg's test p = 0.619, 
Egger's test p = 0.274). The shape of the funnel plot for the homozygote comparison appeared to be 
approximately asymmetrical, which was confirmed by the test results (aa vs. bb, Begg's test p = 0.022, 
Egger's test p = 0.004). However, adjusting the model by the trim and fill method did not infiuence the 
conclusions (OR = 1.23, 95% CI = 1.02-1.48) (Figure 3). 

Figure 3. Begg's funnel plot of publication bias test (homozygote comparison), (a) Before 
trim and fill method; and (b) After trim and fill method. 
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2.5. Discussion 

The meta-analysis revealed that eNOS 4b/a poljmiorphisms were significantly associated with CAD 
in three genetic models (dominant model, homozygote comparison, and heterozygote comparison). 
These data suggest that carriers of the a allele of the eNOS 4b/a gene might be predisposed to CAD. 
Heterogeneity was observed in these three models; however, subsequent meta-regression analyses 
revealed that case and total sample size could explain the heterogeneity. Subgroup analysis based on 
source of the control population revealed that population-based subgroups with the aa or ab genotypes 
had a significantly increased risk of CAD. In stratified analyses according to ethnicity, significantly 
increased risk was detected in the African subgroup with the aa or ab genotypes and the Asian subgroup 
with the aa genotype. Positive associations were also identified in studies whose case or control sample 
size was <500 with the aa or ab genotypes. Similarly, significant associations were identified in studies 
with the aa or ab genotypes in a subgroup analysis according to MI. 

In 1998, a Japanese study including 413 subjects reported that the plasma levels of mono-nitrogen 
oxides in the a allele group (31.2 ± 2.00 |jmol/L) were significantly lower than in the b allele group 
(35.5 ± 0.93 ixmol/L) [48]. Reduced plasma concentrations of nitrogen oxides (NO) in carriers of the a 
allele were also identified in a study by Rittig et al. in a German population [49]. Prevailing experimental 
and clinical data suggest that decreased NO bioavailability accelerates the progression of atherosclerosis 
presumably through mechanisms such as platelet activation, vascular smooth muscle proliferation, 
leukocyte adhesion to the endothelium and increased vascular production of reactive oxygen 
species [50,51]. Therefore, it is possible that eNOS 4b/ a gene polymorphisms could contribute to CAD risk 
by decreasing NO production, which subsequently affects the development of CAD. However, because this 
variant is located in the intronic region, it is possible that it is in linkage disequilibrium with other 
functional variants in regulatory regions of the NOS gene. In 2002, Wang et al. demonstrated that the 
eNOS 4a allele coordinated with the T-786C variant in the promoter region to regulate the transcriptional 
efficiency of eNOS in a haplotype-specific fashion [52]. In 2013, Name et al. reported that the eNOS 
4b/a polymorphism affected the bioavailability and activity of NO via linkage disequilibrium with 
T-786C and G894T variants [25]. Consistent with these observations, our meta-analysis revealed that 
individuals carrying the eNOS 4a allele had a higher CAD risk than subjects carrying the 4b allele. 

Our results also suggested that the 4a allele might be a risk factor for CAD among population-based 
but not hospital-based subgroups. One explanation for this is that there were differences in the frequency 
of exposure to risk factors between hospital- and population-based studies. The use of hospital-based 
population might limit the general applicability of the data. Furthermore, the sample size and number of 
studies in the hospital-based group might not be sufficient to evaluate potential associations. 

In the subgroup analysis according to ethnicity, an increased risk in 4a carriers was identified 
among Afi-ican and Asian populations but not Caucasians. The results could be explained by varying 
susceptibility of different ethnicities to the eNOS 4b/ a poljmiorphism among different ethnicities. 
Environmental characteristics might also contribute to this discrepancy. 

Two previous meta-analyses assessing eNOS 4b/a and CAD risk have been performed. In 2004, 
one study reported an association between eNOS 4a/a and ischemic heart disease compared with 
eNOS 4b/b [10]. Although IHD and CAD have overlapping pathogeneses, they are different disease. 
A second meta-analysis performed in 2010 did not reveal any association between this eNOS genetic 
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variant and CAD risk [11]; therefore, our results contradict this study. This could be explained by several 
reasons. First, different inclusion and exclusion criteria were used. We excluded four studies in 
which the control population deviated from Hardy-Weinberg equilibrium (HWE) [53-56]. We also 
incorporated novel data published between 2010 and October 2013 in our study. Therefore, the number 
and content of the studies included by Li et al. and in the current meta-analysis differed. Second, we 
performed overall analyses in four genetic models; in contrast Li et al. did not perform homozygote 
comparison and heterozygote comparisons. We also performed subgroup analysis based on ethnicity, the 
source of the control population, and case sample size to identify all possible associations between 4b/a 
polymorphisms and CAD; these factors were not investigated previously. 

The present study has some limitations. First, heterogeneity was detected in three genetic models, 
suggesting that results should be interpreted with caution. Second, most of the included studies were 
case-control studies, where the cases were survivors of cardiovascular events; those who did not survive 
were not enrolled. Finally, the sample size in the studies was inadequate. The insufficient sample size 
might increase the probability of false positives or false negatives. 

3. Experimental Section 

3.1. Literature Search 

We systematically searched PubMed using the following medical subject headings (MeSH) or 
search terms: ("CHD" OR "atherosclerosis" OR "coronary atherosclerosis" OR "aortic atherosclerosis" 
OR "coronary heart disease" OR "coronary artery disease" OR "coronary disease" OR "myocardial 
infarction" OR "MI" OR "acute coronary syndrome") and ("polymorphism*" OR "variation") and 
("eNOS" OR "endothelial nitric oxide synthase"). The search covered all English language 
publications within a range of published years from 1996 to 2013 (the last search was performed on 
31 October 2013). The bibliographies of previous meta-analyses and published reviews were also 
checked for additional relevant publications. 

3.2. Inclusion and Exclusion Criteria 

The inclusion criteria for studies were as follows: (1) independent case-confrol, cohort or 
cross-sectional studies; (2) studies that evaluated the association between eNOS infron 4b/a variable 
number of tandem repeats (VNTR) polymorphisms and CAD susceptibility; (3) studies that provided 
complete data regarding genotype number and allele frequencies; (4) studies in which CAD was 
diagnosed based on coronary angiography as well as clinical criteria with clearly reported details; 
(5) studies that also included confrols demonstrated to lack CAD either by coronary arteriography or 
by clinical criteria; and (6) studies that were published in English. The exclusion criteria were as 
follows: (1) insufficient data regarding genotypes and allele frequencies; (2) unconfirmed diagnosis 
of CAD; (3) the genotype disfribution of the control population did not conform to HWE; and 
(4) overlapping publications. 
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3.3. Data Extraction 

Two researchers (Y.Y. and K.D.) extracted information including the first author, year of 
pubhcation, country, ethnicity of the study population, source of controls, genotype number in cases 
and controls, HWE in controls and study outcome. Data were extracted separately and compared to 

reach a consensus. Different ethnicities were categorized as Asians and non-Asians (Caucasians, 
Africans, or others). Different ethnic descents in one study were considered to be independent studies 
when subanalyses based on ethnicity were performed. 

3.4. Statistical Analysis 

The association between eNOS 4b/a gene poljmiorphisms and CAD was compared using the 
odds ratio (OR) corresponding to a 95% confidence interval (95% CI). We used four different genetic 
models in our analysis: homozygote comparison {aa vs. bb), heterozygote comparison {ab vs. bb), 
dominant {aa + ab vs. bb), and recessive models {aa vs. ab + bb). 

Heterogeneity between studies was assessed using p values for the Q statistic and / values. 
Significance was set at the p <0.\ level; p <0.\ for the Q test indicated there was heterogeneity across 
the studies and /? > 0.1 for the Q test indicated a lack of heterogeneity among studies. 
The random-effects model was used to estimate the pooled OR according to the DerSimonian and 
Laird method [57]. / was calculated to describe the percentage of variation caused by the 
heterogeneity: 0%-25%, no heterogeneity; 25%-50%, moderate heterogeneity; 50%-75%, large 
heterogeneity; and 75%-100%, extreme heterogeneity). 

To determine the source of heterogeneity across studies, logistic meta-regression analyses were 
performed. Publication year, population ethnicity, source of the control group, outcome, total sample 
size, case sample size, control sample size, and the ratio of case size to control size were examined to 
evaluate the extent to which different variables explained any observed heterogeneity. Stratified 
analyses were also performed for ethnicity, source of controls, and case sample size. Sensitivity 
analyses were conducted by removing each study individually and consecutively from the combined 
data set. Re-analysis of the remaining studies, allowed the robustness of the main findings to be tested. 
In addition, a funnel plot was used to estimate the potential publication bias; Begg's and Egger's tests 
were then used to examine funnel plot asjmimetry with significance set at < 0.05, and the trim and 
fill method was performed when bias was detected. STATA 12.0 software was used to perform all 
statistical analyses (StataCorp, College Station, TX, USA). 

4. Conclusions 

In conclusion, the results of the current meta-analysis suggested that eNOS 4b/a poljmiorphisms 
could increase the risk of CAD, particularly in African populations and among the population-based 
studies. Further studies with larger sample sizes should be performed to confirm these findings. 
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